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Abstract. A high accuracy photometry algorithm is needed to take full advantage of the potential of the transit 
method for the characterization of exoplanets, especially in deep crowded fields. It has to reduce to the lowest 
possible level the negative influence of systematic effects on the photometric accuracy. It should also be able to cope 
with a high level of crowding and with large scale variations of the spatial resolution from one image to another. A 
recent deconvolution-based photometry algorithm fulfills all these requirements, and it also increases the resolution 
of astronomical images, which is an important advantage for the detection of blends and the discrimination of 
false positives in transit photometry. We made some changes to this algorithm in order to optimize it for transit 
photometry and used it to reduce NTT/SUSI2 observations of two transits of OGLE-TR-113b. This reduction 
has led to two very high precision transit light curves with a low level of systematic residuals, used together with 
former photometric and spectroscopic measurements to derive new stellar and planetary parameters in excellent 
agreement with previous ones, but significantly more precise. 

Key words, planetary systems - stars: individual: OGLE-TR-113 - techniques: image processing - techniques: 
photometric 

1. Introduction reduction method, even with large telescopes, because of 

the high level of crowding present in most of the deep 
Among the -200 exoplanets known so far, only the 10 ones fields of view in the Ga i act i c plane. Nevertheless, the 
C3 ■ transiting their parent star have measured masses and accura t e photometric monitoring of their transits is im- 
radii, thanks to the complementarity of the radial-velocity portant to better constrain the mass-radius relationship 
and transit methods. Among them, 5 were detected by the of close _ in giant planets, and thus the processes of planet 
OGLE-I II . plan etary transit survey (|Udalski et al 2002a | f ormation5 migration and evaporation. Besides, high 
El El |2003J: OGLE-TR-lOb ^on acki et al. 200 5j ), accuracy transit observations may allow the detection 
OGLE-TR-56b ( Kona ckj_et_aL_2 003 , Bou chy et al. 2005| ), of other planetS; ev en terrestrial ones in the best cases, 
OGLE-TR-11 lb ( l E ont et al 20 0^, O GLE-TR-113b by the measurements of the dynamically induced varia- 
( Bouchv e FaT. 2004| IKonacki et al. 2004 J and OGLE- tions of the period of the transit lIMIrllda-Escude 20021 
TR-132b HBouchy et al. 20041 Compared to the other |Agol et al. 20051 |Holman & Murray 20051 . 



transiting exoplanets, they orbit much fainter stars, 
leading to a lower amount of information available from 

their observation. Furthermore, obtaining high accuracy An image deconvolution algorithm 

photometry for these stars is difficult with a classical jMagain et al. 1998| ) has recently been adapted 

to the photometric analysis of crowded fields 
Send offprint requests to: M. Gillon JMagain et al. 20061), even when the level of crowd- 
email : michael.gillon@obs.unige.ch m S * s so high that no isolated star can be used to obtain 
* Based on observations collected with the SUSI2 imager at the PSF (Point Spread Function). We made some mod- 
the NTT telescope (La Silla Observatory, ESO, Chile) in the ifications to this algorithm to optimize it for follow-up 
programme 075.C-0462A. transit photometry, with a main goal in mind: to obtain 
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the highest possible level of photometric accuracy, even 
for faint stars located in deep crowded holds. 

This new method was tested on new photometric ob- 
servations of two OGLE-TR-113b transits obtained with 
the NTT/SUSI2 instrument. This planet was the second 
one confirmed from the list of planetary candidates of the 
OGLE-III survey. It orbits around a faint K dwarf star (J 
= 14.42) in the constellation of Carina. Due to the small 
radius of the parent star (R ~ 0.8 Rq), the transit dip 
in the OGLE-III light curves is the largest one among the 
planets detected by this survey (~ 3 %). As OGLE-TR- 
113 lies in a field of view with a high level of crowding, this 
case is ideal to validate the potential of our new method. 

Sect. 2 presents the observational data. Sect. 3 sum- 
marizes the main characteristics of the deconvolution al- 
gorithm and describes the improvements we brought to 
optimize it for follow-up transit photometry. In Sect. 4, 
our results are presented and new parameters are derived 
for the planet OGLE-TR-113b. Finally, Sect. 5 gives our 
conclusions. 

2. Observations 

The observations were obtained on April 3rd and 13th, 
2005 with the SUSI2 camera on the ESO NTT (pro- 
gramme 075.C-0462A). In total, 235 exposures were ac- 
quired during the hrst night, and 357 exposures during 
the second night, in a 5.4' x 5.4' field of view. The ex- 
posure time was 32 s, the read-out time was 23 s, and 
the i?#813 filter was used for all observations. We used 
SUSI2 with a 2 x 2 pixel binning in order to get at the 
same time a good spatial and a good temporal sampling. 
The binned pixel size is 0.16". The measured seeing varies 
between 0.85" and 1.38" for the first night and between 
1.29" and 1.81" for the second night. Transparency was 
high and stable for both nights. The airmass of the field 
decreases from 1.35 to 1.18 then grows to 1.20 during the 
first sequence, and decreases from 1.23 to 1.18 then grows 
to 1.54 during the second sequence. 

The frames were debiassed and flatfielded with the 
standard ESO pipeline. 

In addition to these new data, we used VLT-FLAMES 
radial velocity measurements QBouchy et al. 2004| ), 
stellar parameters derived from VLT-UVES spec- 
tra (Santo s et al. 200"6jl and OGLE-III ephemeris 
IjKonacki et al. 2004JI to constrain the physical and 
orbital parameters of OGLE-TR-113b. 

3. Photometric reduction method 

3.1. MCS deconvolution algorithm 

The MCS deconvolution algorithm | |Magain et al. 1998| 
hereafter Ml) is an image processing method specially 
adapted to astronomical images containing point sources, 
which allows to achieve (1) an increase of the angular res- 
olution, (2) an accurate determination of the positions 



(astrometry) and the intensities (photometry) of the ob- 
jects lying in the image. One of its main characteristics 
is to perform a partial deconvolution, in order to obtain 
a final image in agreement with the sampling theorem 
l|Shannon 19491 IPress et al. 1989|> . This partial deconvo- 
lution is done by using, instead of the total PSF, a partial 
PSF which is a convolution kernel connecting the decon- 
volved image to the original one. 

In Ml, the determination of the partial PSF was 
not thoroughly addressed. When an image contains suf- 
ficiently isolated point sources, their shape can be used to 
determine an accurate PSF. However, this simple PSF 
determination is rarely possible in crowded fields, which 
generally contain no star sufficiently isolated for this pur- 
pose. 

Magain et al. (2006, hereafter M2) have thus developed 
a version of the algorithm allowing to simultaneously per- 
form a deconvolution and determine an accurate PSF in 
fields containing exclusively point sources, even if no iso- 
lated star can be found. It relies on the minimization of 
the following merit function: 

- 1 

S = J2^( d i-l s *fk) 2 + XH ( s ) (!) 

where * stands for the convolution operator, N is the num- 
ber of pixels within the image, di and <7j are the measured 
intensity and standard deviation in pixel i, Sj is the un- 
known value of the partial PSF and fi is the intensity of 
the deconvolved image in pixel i. H(s) is a smoothing con- 
straint on the PSF which is introduced to regularize the 
solution and A is a Lagrange parameter. This algorithm 
performs an optimal PSF determination, in the sense that 
it uses for this purpose the whole information available in 
the image. It relies on the assumption that the PSF is 
constant over the image. To extend the validity of this 
assumption, one can treat relatively small sub-images if 
PSF variations are suspected. The PSF determination is 
decomposed in several steps and is optimized in order to 
avoid including faint blending stars in the PSF wings, al- 
lowing their detection after inspection of the deconvolved 
image and the residuals map (see M2 for more details). 
Taking into account the blending stars which are unde- 
tectable in the original image results in a better accuracy 
on the PSF, and thus on the astrometry and photometry. 
The deconvolved light distribution / may be written: 

M 

f( x ) = ^2a k r(x - c fc ) (2) 

k=l 

where M is the number of points sources in the image, 
r(x) is the final PSF (fixed) while a k and are free 
parameters corresponding to the intensity and position of 
point source number k. Note that the right-hand side of 
(0) represents only point sources, thus the sky background 
is supposed to be removed beforehand, and it is assumed 
that the data do not contain any extended source. 
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3.2. Optimization of the algorithm for transit 
photometry 

Increase of the processing speed If we consider an image 
with N pixels, containing M point sources, we are left 
with the problem of determining N + 3M parameters, 
i.e. N pixel values of the partial PSF and 3 parame- 
ters for each point source (one intensity and two coor- 
dinates). In follow-up transit photometry, we have gen- 
erally to analyze several hundreds of images for a single 
transit. Furthermore, we are not allowed to analyze only a 
small fraction of the image around the target star. Indeed, 
several systematic noise sources exist, mainly due to at- 
mospheric effects, and the correction of the light curve of 
the analyzed star by the mean light curve of several com- 
parison stars is needed to tend towards a photon noise 
limited photometry. This implies that the S/N (signal-to- 
noise ratio) of the comparison light curve must be signif- 
icantly higher than the S/N of the target star. Thus, we 
have to analyze a field of view large enough to contain 
many reference stars (but yet smaller that the coherence 
surface of the systematics). In practice, we thus need to 
process several hundreds of images containing dozens or 
even hundreds of point sources each. As the deconvolution 
of such an image with the algorithm presented in M2 can 
last up to one day for a very crowded field with an up- 
to-date personal computer, we have to increase drastically 
the processing speed. 

To reach this goal, we use three bits of prior knowl- 
edge. 

- First, we know that our hundreds of images correspond 
to the same field of view and, thus, contain the same ob- 
jects. 

- Secondly, we can assume that the stellar positions do not 
change during the observing run. As we are not interested 
in the astrometry of the stars but only in their photom- 
etry, we can use the best seeing image or a combination 
of the best-quality images as a reference frame and ana- 
lyze it with the standard algorithm in order to obtain the 
astrometry, which is kept fixed during the rest of the anal- 
ysis. All we still need to know to obtain the positions of the 
stars in each image is the amount by which this image is 
translated with respect to the reference image. This trans- 
lation is simply determined by a cross-correlation of the 
images. We neglect image stretch as in practice we treat 
several relatively small sub-images to improve the validity 
of a constant PSF assumption. The result is that, for each 
point source, we are left with only one free parameter (its 
intensity) instead of three. 

- The third prior knowledge is that the relative intensities 
of most point sources do not change much from one image 
to another. We can thus obtain a first approximation of 
the partial PSF by assuming that the relative intensities 
of the point sources are identical to those in the reference 
image, just allowing for a common scaling factor on the 
whole image, which takes into account variations of at- 
mospheric transparency, airmass, exposure time, etc. We 
then use this first partial PSF estimate to obtain a better 



approximation of the intensities. With these better inten- 
sities, we redetermine an improved partial PSF, and so 
on until convergence. 

The most time-consuming task in the standard algo- 
rithm is the iterative determination of the point sources' 
positions and intensities. Here, we already save a large 
amount of computing time by keeping the positions fixed. 
Moreover, when the only unknowns are the point sources' 
intensities, the problem becomes linear in all the param- 
eters. We thus have to solve a set of M linear equa- 
tions, which can be done directly, without any iteration. 
However, as the direct solution of this set of equations is 
quite unstable, we use the Singular Value Decomposition 
method (SVD, IPress et al. 1989|) . which has been found 
to give excellent results. 

The analysis of a set of images is thus divided into two 
parts. In the first one, a reference image is deconvolved 
in order to obtain the astrometry and starting values for 
the point source intensities. Then, the shift of each image 
in the set relative to the reference frame is determined 
by cross-correlation. In the second part, an initial PSF 
is determined for each image, using the fixed astrometry 
(including shift) and photometry. Improved source inten- 
sities are then obtained by solving the linear problem. As 
the accuracy of the partial PSF depends on the accuracy 
of the photometry and vice-versa, the process is repeated 
several times until convergence. In practice, the conver- 
gence is reached after a maximum of 5 cycles. 

Determination of the sky background A tricky problem 
in crowded field photometry is the determination of the 
sky background. Fitting a rather smooth surface through 
seemingly "empty" areas may lead to seeing-dependent 
systematic errors. A much more robust method consists 
in determining the sky background level so that the shape 
of all point sources remains the same, irrespective of their 
intensities and positions. Indeed, a wrong sky level would 
affect weaker sources much more strongly than brighter 
ones. The fact that our method forces all point sources to 
have the same PSF shape can thus be used to obtain an 
accurate determination of the sky background. 

In practice, this is very simply done by not subtract- 
ing the sky background prior to processing, but rather by 
implementing its determination into the method. In this 
case, the observed light distribution d can be modelled as: 

M 

d(x) = s(x) * ^2 a>kr(x - c fe ) + b(x) (3) 
k=l 

where the sky background is represented by the function 
b(x), chosen to be relatively smooth. A 2-dimension sec- 
ond order polynomial (6 free parameters) was found suit- 
able for images obtained in the optical. 

For the deconvolution of the reference image, we have 
now to minimize the following merit function: 

N 1 

5 = E "7 ( d i - bi- [s * m 2 + XH(s) (4) 
i=i Qi 
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where hi is the sky level in pixel i. 

For the deconvolution of the complete set of images, 
the coefficients of b(x) are determined by adding an extra 
step to the analysis of the whole set of images, in each 
iteration. Using the previous approximation of the partial 
PSF and point source intensities, we determine the poly- 
nomial coefficients of the sky background by simple SVD 
solution of a linear set of equations where all parameters 
are fixed but the coefficients of b(x). The whole process 
(f: sky background, 2: partial PSF, 3: intensities deter- 
mination) is repeated until convergence. 

4. Results 

4.1. Light curve analysis 

The light curves obtained with our deconvolution-based 
photometry algorithm are shown in Fig. 1. The flux vari- 
ations before the second transit are intruiging, but we re- 
marked that they are correlated to the location of a bad 
column of the CCD close to OGLE-TR-f 13 and a bright 
reference star. In fact, the bad column is located on their 
PSFs in the first 101 images, and it moved away drasti- 
cally a few exposures before the transit, so the rest of the 
light curve is reliable. The 101 first points were not used 
in the transit fitting. 

For the first night, the dispersion of the light curve 
of OGLE-TR-113 before the transit is 1.20 mmag, while 
the mean photon noise is 0.95 mmag. For the second 
night, the dispersion of the light curve after the transit is 
1.26 mmag, for the same mean photon noise. The slightly 
higher dispersion for the second night can be explained 
by the increased seeing and the fact that OGLE-TR-113 
has a 0.4 mag brighter visual companion about 3" to the 
South (see Fig. 2). When a star's PSF is blended with 
another one, a part of the noise of the contaminating star 
is added to its own noise, resulting in a decrease of the 
maximal photometric accuracy attainable. This effect is 
of course very dependent on the seeing, and may have 
a large impact on the final harvest of a transit survey 
(Gillon & Mag ain, in prep.fr . As the average seeing was 
higher during the second night, we thus expect a lower 
accuracy for this sequence. Nevertheless, the obtained ac- 
curacies for both nights can be judged as excellent. 

Our method has the advantage to produce higher res- 
olution images which can be used to detect a faint blend- 
ing companion around a star which could not be seen on 
a lower resolution image. As shown in Fig. 2, there is no 
evidence of such faint companions around OGLE-TR-113 
in our results. 



4.2. Transit fitting 

The transit fitting was performed with transit curves 
computed with the procedure of Mandel & Agol ( 2002 ) , 
using quadratic limb-darkening coefficients. The transit 
parameters were obtained in two iterations. A prelim- 
inary solution was first fitted to determine epochs for 
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Fig. 1. Our NTT/SUSI2 light curve for the first (top) and 
the second (bottom) observed transits of OGLE-TR-113b, 
with the best fit transit curve superimposed. For the sec- 
ond transit, the variations of the flux before the transit 
are due to a bad column of the CCD located close to the 
PSF cores of OGLE-TR-113 and a bright reference star 
(open symbols). 



NTT photometric series. A period was then determined 
by comparing these epochs with OGLE-III ephemeris 
i|Konacki et al. 2004|l . allowing a very high accuracy on 
the period due to the large time interval separating the 
two sets of measurements (795 times the orbital period). 
The period was then fixed to this value and the radius ra- 
tio, orbital inclination, transit duration and transit epoch 
were fitted by least squares using the NTT data. The limb- 
darkening coefficients used were u± — 0.55 and ui — 0.18, 
obtained from Claret (2000) for the following stellar pa- 
rameters: effective temperature T e ff — 4750 K, metallicity 
[M/H] = 0.1, surface gravity log g — 4.5 and microtur- 
bulence velocity £ t = 1.0 km/s, based on the parameters 
presented in Santos et al. (2006). 

To obtain realistic uncertainties for the fitted transit 
parameters, it is essential to take into account the corre- 
lated noise present in the light curves, as shown by Pont 
et al. p006, in prep.| ). Although we have attained a very 
good level of stability in our photometry, the residuals are 
not entirely free of covariance at the sub-millimag level. 
We model the covariance of the noise from the residuals 
of the light curve itself. We estimate the amplitude of sys- 
tematic trends in the photometry from the standard devi- 
ation over one residual point, o\ , and from the standard 
deviation of the sliding average of the residuals over 10 
successive points, cio- The amplitude of the white noise 
a w and the red noise a r can then be obtained by resolution 
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Fig. 2. OGLE-TR-113 (marked with a cross) in a 256 pixels x 256 pixels sub-image (0.7 ' x 0.7 ') from the worst 
(left) and best (middle) seeing NTT/SUSI2 image of the run (top = North, left = East). The nearby star just South 
of OGLE-TR-113 is about 0.4 mag brighter. Right: deconvolved image. At this resolution (2 pixels = 0.32"), no trace 
of other faint companions is visible. 



of the following system of 2 equations: 

°i = °l + °r (5) 
2 

10 = To r ' ' 

We obtained oy — 400 /imag for both nights. We as- 
sume that a systematic feature of this amplitude could 
be present in the data over a length similar to the transit 
duration. Therefore, if Xbf is the x 2 01 the best fit, in- 
stead of using Ax 2 = 1 to define the 1-sigma uncertainty 
interval, we use 

2 

A X 2 = 1 + N trtl ^- (7) 

for each individual transit, where N tr ,i is the number of 
points in the transit i. As the residuals between different 
transits are not correlated, combining the data from both 
individual transits gives: 

A X 2 = y/Axl + Axl (8) 

Exploring the parameters space for x 2 = Xbf + ^x 2 , 
we then estimated the uncertainties on our parameters. 

The results are given in Table Q The fit of the final 
light curve on the NTT data is shown in Fig. 3. This 
figure also shows the final transit curve superimposed on 
the OGLE-III data (the limb-darkening coefficients are 
changed, as OGLE-III observations have been obtained 
in the I filter), and the phased NTT/SUSI2 data after 
binning on 2 points superimposed on the best-fit transit 
curve. For this binned light curve, the dispersion before 
the transit is ^800 /imag (first night), and ^850 /xmag 
after the transit (second night). 

4.3. Radius and mass determination 

Combining the constraints from our new transit curves 

—1/3 

(dependence on i?»Af» ) and from the spectroscopic de- 
termination of T e ff, log g and [Fe/H] l|Santos et al. 20 06). 



To,ogle[BJD] 2452325.79823 ± 0.00082 
To,! [BJD] 2453464.61665 ± 0.00010 
r , 2 [BJD] 2453474.64348 ± 0.00017 

Table 2. Ephemeris of OGLE-TR-113b transits. T , O gle 
is from Konacki et al. (2004), while To,x and Tq^ are for 
the two transits analyzed in this work. 



we computed the radius and mass of OGLE-TR-113 as in 
Bouchy et al. (2005), taking the relation between M*, 
and the atmospheric parameters from an interpolation of 
Girardi et al. I|2002|l stellar evolution models. The value 
of the planetary radius was then derived from the radius 
ratio and from Next, we fitted a sinusoidal orbit by 
least-squares to the radial velocity data with the new pe- 
riod and epoch, obtaining the planetary mass from M* 
and the semi-amplitude of the radial velocity orbit. 

Our values for the radius and mass of OGLE-TR-113 
and its planetary companion are given in Table ^ which 
also presents the values obtained by Bouchy et al. (2004) 
and Konacki et al. (2004). Our results are in good agree- 
ment with the previous studies, but the uncertainties on 
the mass and on the radius of the planet are significantly 
lower. In fact, our high photometric accuracy allows to 
reach the regime where the uncertainties on the mass and 
radius of the primary dominate: the use of other stel- 
lar evolution models should introduce parameter changes 
which are of the same order than the error bars. The de- 
velopment of this point is beyond the scope of this paper. 

4.4. Transit timing 

OGLE-III <|Konacki et al. 2004jl and our new NTT tran- 
sits ephemeris are presented in Table |21 Using these 
ephemeris and assuming the absence of long-term transit 
period variations due, e.g., to the oblateness of the star or 
dissipative tidal interactions between the planet and the 
star, we obtained a value for the period in perfect agree- 
ment with the one presented in Konacki et al. (2004), but 
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Fig. 3. The best-fit transit curve is shown together with the phased NTT/SUSI2 data after binning on 2 points (top), 
without binning (middle) and with the phased OGLE-III data (bottom). 



with a much smaller error bar. Our extremely high ac- 
curacy on the orbital period (~ 0.1 s) is mainly due to 
the long delay between OGLE-III and NTT observations. 
Moreover, our accuracy on the epoch of the NTT transits 
is ~ 12 s, while the accuracy on the OGLE-III epoch is ~ 
71 s. 

Short-term transit timing variations (TTV) 
may be induced by the presence of a satel- 
lite IjSartoretti fc Schneider 1999jl or a second 
planet (|Miralda-Escude 20021 |Agol et al. 20051 
|Holman fc Murray 2005| ). If we examine the temporal 



deviation of the second NTT transit off the expectation 
based on our period and the first NTT transit time, we 
obtain a value of ~ 43 s. This TTV has a statistical 
significance of 2.5 sigmas. Nevertheless, we must stay 
cautious because systematics are able to slightly distort 
light curves, all the more so since these transit light 
curves are not complete. This is the case here: we lack 
the flat part after the first transit and, due to the 
bad column, the flat part prior to the second one. A 
way to estimate the likelihood of the hypothesis that 
systematics are responsible for the observed TTV is to 
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A B C 

Inclination angle [deg] 88.8 - 90 85 - 90 88.4 ± 2.2 

Period [days] 1.4324757 ± 0.0000013 1.43250 (adopted) 1.4324758 ± 0.0000046 

Semi-major axis[AU] 0.0229 ± 0.0002 0.0228 ± 0.0006 0.02299 ± 0.00058 

Eccentricity (fixed) 

Primary mass[M Q ] 0.78 ± 0.02 0.77 ± 0.06 0.79 ± 0.06 (adopted) 

Primary radius [R e ] 0.77 ± 0.02 0.765 ± 0.025 0.78 ± 0.06 (adopted) 

Planet mass [Mj] 1.32 ±0.19 1.35 ±0.22 1.08 ±0.28 

Planet radius[flj] 1.09 ± 0.03 l-08-o°ol 1.09 ±0.10 

Planet density[g cm -3 ] 1.3 ± 0.3 1.3 ± 0.3 1.0 ± 0.4 

Table 1. Parameters obtained from this analysis (A) for OGLE-TR-113 and its planetary companion, compared to 
the ones presented in Bouchy et al. (2004) (B) and in Konacki et al. (2004) (C). 



^0,1, ingress [BJD] 
7b, 1, egress [BJD] 

Tb, 2, ingress [BJD] 
Tb, 2, egress [BJD] 



2453464.61669 ± 0.00014 
2453464.61653 ± 0.00014 
2453474.64299 ± 0.00024 
2453474.64400 ± 0.00024 



Table 3. Ephemeris obtained from ingress and egress of 
both transits. 
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Fig. 4. This figure shows the agreement between the tran- 
sit times obtained from the complete light curves (top) and 
using ingress (I) and egress (E) independently (bottom) 
and the predicted transit times based on the determined 
period and OGLE-III ephemeris. 



determine the transit ephemeris for the ingress and egress 
independently. Table and Fig. 4 show the result. For 
both transits, transit times obtained from ingress and 
egress are in good agreement, leading us to reject such 
systematic errors as the explanation for the observed 
TTV. 

The statistical significance of the TTV is rather low, 
as can be seen in Fig. 4. Future observations are needed 



to confirm its existence. We notice nevertheless that our 
timing accuracy would clearly be good enough to allow 
the detection of a third planet or a satellite giving rise 
to a TTV amplitude of 1 minute. Considering the case 
of OGLE-TR-113b and, as the cause of such a TTV, a 
satellite with an orbital distance to the planet equal to 
the Hill radius, we can obtain an estimate of its mass by 
using the formula USartoretti fc Schneider 19 99): 



M, 



St /3M» 



P„ \ M„ 



1/3 



(9) 



where M* , M p and M s are respectively the masses of the 
star, planet and satellite, P p is the orbital period of the 
planet and St the amplitude of the TTV. We obtain M s ~ 
7M e . 

For an exterior pertubing planet, the most interesting 
case would be a pertubing planet in 2:1 mean-motion res- 
onance with OGLE-TR-113b, for which we obtain with 
the following formula ( |Agol et al. 2005| ): 



M, 



p2 



•pi 



pi 



(10) 



a mass ~ lM ffi . 

These computations demonstrate the interest of high 
accuracy photometric follow-up of known transiting ex- 
oplanets and show that the accuracy obtained with our 
reduction method for high quality data would allow the 
detection of very low mass objects. 

5. Conclusions 

The results presented here show that our new photome- 
try algorithm is well suited for follow-up transit photom- 
etry, even in very crowded fields. After analysis of NTT 
SUSI2 observations of two OGLE-TR-113b transits, we 
have obtained two very high accuracy transit light curves 
with a low level of systematic residuals. Combining our 
new photometric data with OGLE-III ephemeris, spectro- 
scopic data and radial velocity measurements, we have 
determined planetary and stellar parameters in excellent 
agreement with the ones presented in Bouchy et al. (2004) 
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and Konacki et al. (2004), but significantly more precise. 
We notice that the sampling in time, the sub-millimag 
photometric accuracy and the systematics residuals level 
of our light curves would be good enough to allow the 
photometric detection of a transiting Hot Neptune, in the 
case of a small star as OGLE-TR-113. 

We have obtained a very precise determination of 
the transit times, and, combining them with OGLE-III 
ephemeris, we could determine the orbital period with a 
very high accuracy. The precisions on the epochs and the 
period would in fact be high enough to allow the detection 
of a second planet or a satellite, for some ranges of orbital 
parameters and masses. Even a terrestrial planet could be 
detected with such a transit timing precision. 
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